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Abstract 

We investigate the recently proposed natural double inflation model in super- 
gravity. Chaotic inflation first takes place by virtue of the Nambu-Goldstone- 
like shift symmetry. During chaotic inflation, an initial value of second infla- 
tion (new inflation) is set, which is adequately far from the local maximum 
of the potential due to the small linear term in the Kahler potential. Then, 
primordial fluctuations within the present horizon scale may be produced 
during both inflations. Primordial fluctuations responsible for anisotropies 
of the cosmic microwave background radiation and the large scale structure 
are produced during chaotic inflation, while fluctuations on smaller scales are 
produced during new inflation. Because of the peculiar nature of new infla- 
tion, they can become as large as 10 -1 -10 -2 , which may lead to the formation 
of primordial black holes. 

PACS number(s): 98.80.Cq,04.65.+e,12.60.Jv 
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I. INTRODUCTION 



Inflation is the most attractive mechanism to generate primordial density fluctuations 
responsible for anisotropies of the cosmic microwave background radiation (CMB) and the 
large scale structure, in addition to solving the flatness and the horizon problem |E]. Re- 
alistic inflation models should be constructed in the context of the supersymmetric theory, 
especially, its local version, supergravity (SUGRA) 0] because supersymmetry (SUSY) guar- 
antees the flatness of the inflaton against radiative corrections and gives a natural solution 
to the hierarchy problem between the inflationary scale and the electroweak scale [||. 

New inflation [|J is very attractive in the context of SUGRA because it takes place 
at a low energy scale and naturally leads to a sufficiently low reheating temperature to 
avoid the overproduction of gravitinos. However, roughly speaking, new inflation has two 
severe problems One is the initial value problem: that is, the inflaton must be fine- 
tuned near the local maximum of the potential for sufficient inflation. The other is the 
flatness (longevity) problem: that is, why the universe lives so long beyond the Planck time. 
Asaka, Kawasaki, and the present author || found that, due to gravitationally suppressed 
interactions with particles in the thermal bath, the inflaton can dynamically go to the local 
maximum of its potential. However, the other problem still exists unless the universe is open 
at the beginning. Izawa, Kawasaki, and Yanagida || considered another type of inflation 
(called preinflation) which takes place before new inflation and drives the inflaton for new 
inflation dynamically toward the local maximum of its potential. If preinflation is chaotic 
inflation 0, the longevity problem is solved too. 

It, however, is believed to be difficult to realize chaotic inflation in SUGRA. This is 
mainly because a scalar potential in minimal SUGRA has an exponential factor with the 
form exp(\<p\ 2 / Mq) so that any scalar field <fi cannot take a value much larger than the reduced 
Planck scale Mq ~ 2.4 x 10 18 GeV. Several supergravity chaotic inflation models were pro- 
posed by use of functional degrees of freedom of the Kahler potential in SUGRA ||||. But, 
there are no symmetry reasons to have such proposed forms. Recently, Kawasaki, Yanagida, 
and the present author [ITJ proposed a natural model of chaotic inflation in SUGRA by 
use of the Nambu-Goldstone-like shift symmetry. Motivated by this model, Yokoyama and 
the present author |TT[ considered chaotic inflation followed by new inflation, where chaotic 
inflation first takes place around the Planck scale to solve the longevity problem and gives 
an adequate initial condition for new inflation. Thus double inflation has been proposed as 
a solution of initial condition problems of some types of inflation |),13|. (See Refs. [12,13| 
for the initial value problem and its solution of hybrid inflation.) 

On the other hand, double inflation is also motivated by observational results. One 
motivation is to reconcile predicted spectra with observations of the large scale structure 
It is known that a standard cold dark matter (CDM) model in a flat universe with a 
nearly scale-invariant spectrum cannot reproduce the observation of the large scale structure. 
Furthermore, the recent observations of anisotropies of the CMB by the BOOMERANG 
experiment |T5] and MAXIMA experiment [TTJ found a relatively low second acoustic peak. 
Another motivation is to produce primordial black holes (PBHs) |IT|-|T9|| . Massive compact 



halo objects (MACHOs) are observed through gravitational microlensing effects [^(J, which 
are a possible candidate of dark matter. Furthermore, PBHs evaporating now may be a 
source of antiproton flux observed by the BESS experiment or responsible for short 
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gamma ray bursts (GRBs) p2| |. Though many double inflations have been considered, 



they are often discussed in a simple toy model with two massive scalar fields. However, 
a natural double inflation model in SUGRA is recently proposed ||23|| , where there are no 



initial condition problems and the model parameters are natural in the 't Hooft sense |24 



In this model, chaotic inflation takes place first of all, during which an initial value of new 
inflation is dynamically set due to the supergravity effects. It can be adequately far from 
the local maximum of the potential due to the small linear term of the inflaton in the 
Kahler potential. Therefore primordial density fluctuations responsible for the observable 
universe can be attributed to both inflations, that is, chaotic inflation produces primordial 
fluctuations on large cosmological scales and new inflation on smaller scales.^] The energy 
scale of new inflation becomes of the same order as the initial value of new inflation so that 
produced density fluctuations may become as large as the order of unity due to the peculiar 
nature of new inflation, which straightforwardly may lead to PBHs formation. 

In this paper we minutely investigate the recently proposed natural double inflation 
model in supergravity, especially, primordial density fluctuations produced during inflation. 
Then, the PBHs formation is discussed. 



II. MODEL AND DYNAMICS 
A. Model 



In this section we briefly review the double inflation model in supergravity proposed 
recently p3fl . We introduce an inflaton chiral superfield 9) and assume that the model, 



especially, Kahler potential K(§, $*,...) is a function of $ + $*, which enables the imaginary 
part of the scalar component of the superfield $ to take a value larger than the gravitational 
scale, and leads to chaotic inflation. Such a functional dependence of K can be attributed 
to the Nambu-Goldstone-like symmetry introduced in Ref. flop . We also introduce a spurion 
superfield H describing the breaking of the shift symmetry and extend the shift symmetry 
as follows, 



$ -> $ + i CM G , 



$ + i CMq 



(1) 



where C is a dimensionless real constant. Below, the reduced Planck scale Mq is set to be 
unity. Under this shift symmetry, the combination S$ 2 is invariant. Inserting the vacuum 
value into the spurion field, (5) = A, softly breaks the above shift symmetry. Here, the 



1 In Ref. (Tl]] the initial value of new inflation is so close to the local maximum of the potential 
for new inflation that the universe enters a self-regenerating stage [25|,26|. Therefore primordial 
fluctuations responsible for the observable universe are produced only during new inflation. Fur- 
thermore, even if chaotic inflation proposed in Ref. [^] is adopted as preinflation in Ref. Q, the 
same situation occurs, that is, second inflation becomes eternal inflation because the superpotential 
in Ref. [p]] vanishes during chaotic inflation. 
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parameter A is fixed with a value much smaller than unity representing the magnitude of 
breaking of the shift symmetry ([!]). 

We further assume that in addition to the shift symmetry, the superpotential is invariant 
under the U(1)r symmetry because it prohibits a constant term in the superpotential. The 
above Kahler potential is invariant only if the .R-charge of $ is zero. Then, we are compelled 
to introduce another supermultiplet X(x, 6) with its i?-charge equal to two, which allows 
the linear term X in the superpotential. As shown in Ref. []23| ], for successful inflation, the 
absolute magnitude of the coefficient of the linear term X must be at most of the order of 
|A|, which is much smaller than unity. Therefore in order to suppress the linear term of X 
in the superpotential, we introduce the Z 2 symmetry and a spurion field II with odd charge 
under the Z 2 symmetry and zero -R-charge. The vacuum value (II) = v softly breaks the 
Z 2 symmetry and suppress the linear term of X. Then, the general superpotential invariant 
under the shift, U(l)# and Z 2 symmetries is given by 

W = a x UX {l + a 2 (H$ 2 ) 2 + •••}- a 3 X {s$ 2 + « 4 (H$ 2 ) 3 + •••}, (2) 

where we have assumed the i?-charge of H vanish and $, X, and H fields are odd under the 
Z 2 symmetry.0(See table I in which charges for superfields are shown.) Here, a, are complex 
constants of the order of unity. 

After inserting vacuum values of spurion fields 5 and EE, the superpotential is given by 

W = vX {l + a 2 (A$ 2 ) 2 + ■■■}- X {A$ 2 + a 4 (A$ 2 ) 3 + ••-}. (3) 

Here, the complex constants a± and a% are renormalized into v and A. Though the above su- 
perpotential is not invariant under the shift and the Z 2 symmetries, the model is completely 
natural in the 't Hooft's sense [E3J] because we have enhanced symmetries in the limit A and 
v — > 0. As long as [$[ <C |A| -1 /^ higher order terms with a« of the order of unity become 
irrelevant for the dynamics. Therefore we can safely omit them in the following discussion. 
After all, we use, in the following analysis, the superpotential, 

W ~ vX - AX$ 2 (4) 
= vX(l-g<S> 2 ), (5) 

with g = X/v. Though, generally speaking, only a constant can become real by the use of 
the phase rotation of the X field, below we set both constants v and A (g) to be real for 
simplicity]^] 

The Kahler potential neglecting a constant term and higher order terms is given by 



2 The oddness of the spurion field 3 under the Z 2 symmetry implies that it breaks both the shift 
symmetry and the Z 2 symmetry at once. So, we expect that the magnitudes of the breaking of both 
the Z 2 and the shift symmetries are of the same order. We hope that the yet unknown mechanism 
simultaneously gives the spurion field S and II the vacuum values and such a mechanism be realized, 
for example, in the superstring theory. 

3 The dynamics of the general case is discussed in Ref. [23]. 
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K = v 2 {$ + $*) + -($ + $*) 2 + XX*. 
2 



(6) 



Here i> 2 ~ t> is a real constant representing the breaking effect of the Z% symmetry. Here and 
hereafter, we use the same characters for scalar with those for corresponding supermultiplets. 

B. Dynamics 

Now that the Kahler potential and the superpotential are specified, the Lagrangian 
density L($,X) for the scalar fields $ and X is given by 



L($, X) = d^d**®* + d^Xd^X* - V(<S>, X). 



(7) 



The scalar potential V of the chiral superfields X(x, 9) and <3>(x, 9) in supergravity is given 
by 



V = v 2 e K 



\x\ 2 + \xn + \x\ 



-2c/$ + (v 2 + $ + $*)(! -g& 



Now, we decompose the scalar field $ into real and imaginary components, 



Then, the Lagrangian density L(ip,x,X) is given by 

L(<p, X , X) = \d^d^ + U^X + dpX&X* - V(tp, X , X), 
with the potential V(ip, x, X) given by 



(8) 



(9) 



(10) 



V2 



V((f,X,X) = ve 2 exp <{ ( y? + 



\X\ 



x 



{ [l " 9 (<p 2 " X 2 ) + \ 9 V + X 2 ) 2 



;i-|x| 2 + |x| 4 ) 



2<?V + x 2 ) 

-(v 2 + V2 V ) {V2 gtp [ 2 - g{ V 2 - X 2 )}- ^ g 2 VX 2 } 
+ (V2 + V2^) 2 [l - g( V 2 ~X 2 ) + \ 9 2 (V 2 + X 2 ) 2 } } } 



(11) 



Because of the exponential factor, ip and X rapidly goes down to 0(1). On the other 
hand, x can take a value much larger than unity without costing exponentially large potential 
energy. Then the scalar potential is approximated as 



V ~ A 2 + 2 X 2 \X\ 2 



(12) 
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with A = gv. Thus the term proportional to x 4 becomes dominant and chaotic inflation can 
take place. Then, using the slow-roll approximation, the e-fold number N during chaotic 
inflation is given by 

~ X~ 

N ~ — . (13) 
The effective mass squared of ip, m 2 ,, during chaotic inflation becomes 

ml * ^ * 6H 2 » H 2 * (14) 

where H is the hubble parameter at that time. Therefore (p oscillates rapidly around the 
minimum (p m i n so that its amplitude damps in proportion to a~ 3//2 with a being the scale 
factor. Here, the potential minimum for (p, (p m i n , during chaotic inflation is given by 



/V2. (15) 



-v 2 

Thus the initial value of the inflaton (p of second inflation (new inflation) is set dynamically 
during chaotic inflation. 

On the other hand, the mass squared of X, m 2 x , is dominated by 

24 

m\ ~ 2A\ 2 ~ —H\ (16) 

which is much smaller than the hubble parameter squared until x 2 ~ 24 so that X also 
slow-rolls. In order to analyze the dynamics of X, we set X to be real and positive making 
use of the freedom of the phase choice. In this regime classical equations of motion for X 
and x are given by 

3HX~ -m 2 x X, (17) 
3if X ~-AY, (18) 



which leads to 



X 




(19) 



where Xi and x% are the initial values of X and x fields. This relation holds actually if and 
only if quantum fluctuations are unimportant for both x an d X. First of all, for Xi the 
comparison of the magnitude of quantum fluctuations and that of the classical evolution 
during one hubble time shows that quantum fluctuations become dominant if x ~ A -1 / 6 , 
when the universe enters the self- reproduction stage of eternal inflation p5| , |26| . So, we 



consider only the regime with x ^ A -1 / 6 , where the classical equation of motion Eq. ([18] 
is valid. Next, we estimate the amplitude of quantum fluctuations for X, Using the Fokker- 
Planck equation for the statistical distribution function of X based on the stochastic inflation 
method of Starobinsky |27| , the root-mean-square (rms) of quantum fluctuations for X, 
((AX) 2 ), is given by (TJ 
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(20) 



On the other hand, using Xi = 0(1) and Xi ~ A" 1 / 6 , the classical value of X becomes 
X ~ A 1 / 3 ^ 2 . Thus, since AC 1, the amplitude of X becomes much smaller than unity by 
the time x — v^24, when the effective mass squared m 2 x is comparable with H 2 . Thereafter, 
X rapidly oscillates around the origin and its amplitude damps in proportion to a _3//2 even 
more. Thus our approximation that both tp and X are much smaller than unity is consistent 
throughout the chaotic inflation regime. 

As x becomes of order of unity, either the constant term v 2 or the term with v 2 gx 2 
becomes dominant. In the former case, small hybrid inflation takes place, which is followed 
by new inflation. Hence there is no break between chaotic and new inflation. On the other 
hand, in the latter case, x rapidly oscillates around the origin until new inflation starts 
so that there is a break between them, though the scale factor grows twice at most. A 
little numerical calculation shows that if g > 1.1, we have a break between chaotic and new 
inflation. 

Next let us investigate when new inflation starts. The potential with X ~ is approxi- 
mated as 



2 

V(<p, Xi x — o) — v e ~ 2 ex P V 9 + 



x 



J\ 2 , 2 1 ,9 2,9 2 

+ X [9+^ + yX 



V2J 

(21) 



The global minima are given by y? 2 = 2/g and x — 0- The mass squared for ip, m 2 , reads 



m 2 ~-(g -l) + (g + \g 2 )x 2 - (22) 



Thus new inflation begins when x — Xcrit given by 



x - t = -gfirl- (23) 

Once new inflation begins, x rapidly goes to zero because the effective mass squared 
becomes m 2 ~ QgH 2 > QH 2 . Then, for x — and X <tC 1, the potential is given by 

V(<p, x ^ 0, X « 1) ~ v 2 {l - (g - l)$ 2 + 2(g - 1) 2 ^ 2 \X\ 2 + •••}, (24) 

where (p = <p — <p ma x and <p m ax = v 2/[V2(g — !)]■ Thus if g > 1 (A > v), new inflation takes 



place and ip rolls down slowly toward the vacuum expectation value r] = y2/g. 

Before new inflation starts, (p stays at (p m i n , which is different from p) max - Then, the 
initial value of <p, <p>i for new inflation is given by 

V2 9 , ,x 

V ' = ~2~l <25) 
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On the other hand, the amplitude of quantum fluctuations of <p is estimated as btp, 
v/(27ry / 3)- Using the fact that g > 1, we find that quantum fluctuations do not dominate 
the dynamics unless v 2 ^ v. 

The total e-folding number N new during new inflation is given by 



2(9 



V2 



fj 



V2 g 



(26) 



Then the total e-folding number N is given by N = N + N new . We set iV C oBE = 60 for 
simplicity, when the physical wavenumber of the mode (&cobe) corresponding to the Cosmic 
Background Explorer (COBE) scale exits the horizon, that is, /ccobe/ (0(^=60)^(^=60)) — 1- 

In case N new > 60, primordial density fluctuations responsible for the observable universe 
are produced only during new inflation. Otherwise, chaotic inflation produces primordial 
fluctuations on large cosmological scales and new inflation on smaller scales. In this paper 
we consider only the latter case. 

After new inflation, <p oscillates around the global minimum r\ and the universe is domi- 
nated by a coherent scalar-field oscillation of a = ip — rj. Expanding the exponential factor 



e v ^ W = e v \l + 2 V a + ■••), (27) 

we find that a has gravitationally suppressed linear interactions with all scalar and spinor 
fields including minimal supersymmetric standard model (MSSM) particles. For example, 
let us consider the Yukawa superpotential W = i/iDiHSi in MSSM, where Di and Si are 
doublet (singlet) superfields, H is a Higgs superfield, and yi is a Yukawa coupling constant. 
Then the interaction Lagrangian is given by 

C int ~ytriaI%S? + ---, (28) 

which leads to the decay width T given by 

r-E^Vm*. (29) 



Here m a ~ 2^ r g^ev^ 9R v is the mass of a. Thus the reheating temperature Tr is given by 

T R ~0.1y V mf } (30) 

where y = vt- Taking y ~ 1, the reheating temperature Tr is given by 

T R ~ v 3 / 2 < A 3/2 . (31) 

As shown later, the upper bound of A is given by A < 1.2 x 10~ 6 . Hence the reheating 
temperature Tr is constrained as 

T R < 10~ 9 ~ 10 9 GeV, (32) 

which is low enough to avoid the overproduction of gravitinos in a wide range of the gravitino 
mass p8 . 
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III. DENSITY FLUCTUATIONS AND PBHS FORMATION 



A. Density fluctuations 

In this section we investigate primordial density fluctuations produced by this double in- 
flation model. First of all we consider density fluctuations produced during chaotic inflation. 
As shown in the previous section, there are two effectively massless fields, \ and X, during 
chaotic inflation. Using Eq. ( |i~2"D and adequate approximations, the metric perturbation in 
the longitudinal gauge $^ can be estimated as [|2U 



H X 2 
1>a = --pC, - 16^C 2 , 

£7, = H S f, 

C 2 = H k- «) A (33) 

V x x J x 

where the dot represents time derivative, the term with C\ corresponds to the growing 
adiabatic mode, and the term with C2 the nondecaying isocurvature mode. You should 
notice that only \ contributes to growing adiabatic fluctuations. Then, with the fact that 
I < 1, the amplitude of curvature perturbation $^ on the comoving horizon scale at 
X = Xfi is given by the standard one-field formula and reads 

where / = 3/5 (2/3) in the matter (radiation) domination. If N new < 60, the comoving 
scale corresponding to the COBE scale exits the horizon during chaotic inflation. Defining 
-^cobe as the e-folding number during chaotic inflation, corresponding to the COBE scale, 
the COBE normalization requires $a(^cobe) — 3 x 1(U 5 [00 . Then the scale A is given by 



A = 4.2 x 10-^ obi . (35) 



The spectral index n s is given by 



n s ~ 1 - ~ . (36) 

jVcobe 



Since the COBE data shows n s = 1.0 ±0.2 §0|, Y CO be > 15, which leads to A < 1.2 x 10~ 6 . 

Next let us discuss density fluctuations produced during new inflation. In this case, both 
tp and X are effectively massless fields. As with the case of chaotic inflation, the metric 



perturbation in the longitudinal gauge $^ can be estimated as |29] 



C = H (| - f ) 2( 9 - l),p\ (37) 
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In this case, also, only (p contributes to growing adiabatic fluctuations. Then, with the fact 
that X < 1, the amplitude of curvature perturbation $^ on the comoving horizon scale at 
(p = (fN is given by 

$ A (N) ~ — 1= -, — V -^~. (38) 

The spectral index n s of the density fluctuations is given by 

n s ~ 1-4(0-1). (39) 

Here we relate the e-folding number N < N new with the comoving wave number k. 
Notice that the hubble parameterjduring new inflation is much smaller than that during the 
early stage of chaotic inflation [N = O(10)]. Then, the e-folding number N < N new when 
the comoving wave number k = fccoBE e^ m ~ N ' (Nqobe = 60) exits the horizon during new 
inflation is determined by 

kcOBE6 1. (40) 



a (N=m) e( 60 "^) H N 

Using H N ~ H( N=e0 ) ~ A%~ /(2v3), and Eq. flT3]), the correspondence is given 

by 



N = N' + In 



9X 2 ~ 

■^COBE - 



= N' -InAgNcoBE- (41) 
The deviation from the standard correspondence (N = N') is not negligible for A^ CO be = 

o(io).g 

You should also notice that (p ~ v 2 at the beginning of new inflation. Then, since v ~ V2, 
the amplitude of curvature perturbation $^4 can become as large as the order unity, which 
may lead to PBHs formation. 

Before discussing PBHs formation, let us comment on the case with a break between 
chaotic and new inflation. As shown before, we have such a break if g > 1.1. In this case 
some of the comoving wave numbers which exit the horizon during chaotic inflation reenter 
the horizon and again exit during new inflation. Therefore, for such modes, we need to 
compare the amplitude of quantum fluctuations induced during chaotic and new inflation. 
Following the procedure as done in Ref. we can easily show that fluctuations induced 



during chaotic inflation are a little less than newly induced fluctuations during new inflation. 
Thus we conclude that the fluctuations of <p induced in chaotic inflation can be neglected 
when we estimate density fluctuations produced during new inflation. 



4 The hubble parameter during chaotic inflation changes significantly (H oc N). Hence, the de- 
viation from the standard correspondence may be also significant during chaotic inflation. The 
e-folding number N < Nqobe when the comoving wave number k = kcoBE^ 1 ^ 00 ^" 1 ^ exits the 
horizon during chaotic inflation is given by N' = N + ln(iVcoBE/A0- 
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B. Primordial black holes formation 



PBHs have been paid renewed attention to because they may explain the existence of 
massive compact halo objects (MACHOs) |2(J and become a part of cold dark matter. 



Furthermore, PBHs are responsible for antiproton fluxes observed by the BESS experiments 



[ 21] or short gamma ray bursts p2[ . 

Carr and Hawking first discussed PBHs formation and showed that in the radiation 
dominated universe, a black hole is formed soon after the perturbed region reenters the 



horizon if the amplitude of density fluctuations 5 lies in the range 1/3 < 5 < 1 ||31|| . Then, 
the mass of produced PBHs Mbh is roughly given by the horizon mass, 



<-BH 



y/p J VGeV 



O.O66M — - , (42) 



where p and T are the total energy density and the temperature of the universe at formation. 
The horizon scale at formation is related to the present cosmological scale L by 

~ a(T) 1 h 

~ 6.4 x l(T 8 Mpc (^gy) , (43) 

with To ~ 2.7 K the present temperature of the universe. The corresponding comoving wave 
number k = 27r/L is given by 

k ~ 1.0 x 10 8 Mpc" 1 (^y) • (44) 

Assuming Gaussian fluctuations, the mass fraction of produced PBHs (/3 = Pbh/p) is 
given by 

"" (M)eXP (-l^(M))' (45) 

where a(M) is the root mean square of mass variance evaluated at horizon crossing. Bullock 
and Primack |32] pointed out that the standard Gaussian assumptions may be inadequate 
because PBHs are formed at high density peaks so that the linear theory may be invali- 
dated and non-Gaussianity may affect the abundance of PBHs significantly (see also [|33]| ). 
However, for example, in the model adopted in [Q, the root mean square of mass vari- 



ance calculated in the standard Gaussian theory differs at most by the factor 1.5 from that 
calculated taking into account the non-Gaussian effects. The aim in this paper is just to 
demonstrate that our double inflation model straightforwardly leads to PBHs formation, 
and the concrete values of the parameters should not be taken seriously. Therefore, for our 
purpose, it is sufficient to assume that fluctuations are Gaussian distributed. The analysis 



of non-Gaussianity in a similar model as ours is done in Ref. [IS . 
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Using the mass fraction (3 at formation, the ratio of the present energy density pbh{M) 
of PBHs with the mass M and the entropy density is given by 

~ ~/?(M)T, (46) 

S 4r 

which yields the normalized energy density, 

T 



n BH h- 2 ~ 2.1 x 10 8 /3 ( 
~ 5.4 x 10 7 /3 ( 



GeV 



MY 1 ' 2 



^MiM- (47) 
Then, MACHO PBHs with mass ~ O.1M are produced at the temperature given by 

T ~ 0.81 GeV, (48) 

which corresponds to 

L ~ 7.9 x 10" 8 Mpc, 

fc~8.1x 10 7 Mpc" 1 . (49) 

As easily seen from Eq. fl38|), the fluctuations with the largest amplitude are produced at the 
onset of new inflation, which we identify with the formation time of PBHs. As shown later, 
the spectrum is so steep that the formation of the PBHs with smaller masses is suppressed 
strongly. From Eq. (|||) we obtain N' new ~ 37, which corresponds to 

iV new ~37-ln92(7. (50) 

On the other hand, the present energy density Qsnh' 2 ~ 0.25 is explained if the mass 
fraction is given by 

/5~1.5xl0" 9 , (51) 

which implies the mass variance a ~ 0.056 under the Gaussian approximation, corresponding 
to 

$4 ~ 0.04. (52) 

Taking into account Eqs. (0), (|52|), and the COBE normalization (j35|), MACHO PBHs are 
produced in this model if we take the parameters given by 

A ~ 1.3 x HT 6 , 
v ~ 1.1 x 10~ 6 , 

v 2 ~ 0.93 x 10~ 6 , (53) 

with g = X/v ~ 1.2. 
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As another example, let us consider PBHs responsible for antiproton fluxes observed by 
the BESS experiments [21] or short gamma ray bursts p2 |. Such PBHs are evaporating 



now, which leads to the initial mass M ~ 3 x 1O~ 19 M . Then, the temperature at formation 
is given by T ~ 4.7 x 10 8 GeV corresponding to L ~ 1.4 X 10~ 16 Mpc and k ~ 4.7 x 10 16 
Mpc -1 . On the other hand, the abundance is given by VtBuh 2 ~ 2 x 10 -9 , which implies 
(3 ~ 2.0 x 10~ 26 , cr ~ 0.032, and $a — 0.002. The PBHs satisfying the above conditions are 
produced if we take the parameters given by 

A ~ 5.8 x 10~ 7 , 
«~3.9x 10~ 7 , 

v 2 ~ 5.2 x 10~ 6 , (54) 

with g = X/v ~ 1.5. Note that in both cases, all parameters are of the same order. Also, 
the temperatures at formation are lower than the reheating temperature so our assumption 
that PBHs are formed in the radiation dominated universe is justified. 



C. Numerical calculations of density fluctuations 

In this section we numerically calculate density fluctuations produced during double 
inflation in order to confirm the analytic results given above and show explicitly that den- 
sity fluctuations which lead to PBHs are realized in our model. Our method of numerical 



calculations is based on Ref. |34| 



We decompose multiscalar fields <f>i(x,t) into the homogeneous mode <j>i(t) and fluctua- 
tions Xi(x, t), 

<f)i{x,t) = 4>i{t) + Xi(x,t) 

= m + J j0T2 ^ k ' x + ^ t ^ lkx ) ■ ^ 

The metric is also expanded around the background metric, 
g^(x, t) = g^(t) + y^x, t) 

= 9,M + / J0r 2 (y^(*, t)e ik - x + yl(k, t)e~ ik x ) . (56) 

Since fluctuations are generated quantum mechanically, we need to treat them as quantum 
Heisenberg operators. Defining creation and annihilation operators satisfying the canonical 



commutation relations [a;(fc), aj(fc')] = SijS(k — k'), fluctuations are expanded over such 
operators 



Xi(k,t) =J2 Sl Ptj( k ^ t )a j (k), 

j 

y^(k,t) =^2h tlVji (k,t)ai(k). (57) 

i 

In particular, the metric perturbation in the longitudinal gauge is expanded as 
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d 3 k 



^ X ^~J (2^ 

$(M) = J2Mk,t)Mk) 



$(k,t)e lkx + &(k,t)e- ikx 



(58) 
(59) 



The equations of motion for the homogeneous mode are given by 

- - dV 

<P t + 3tf & + = 0, 



(60) 



with the hubble constant H given by 



H 2 



3M 2 , 



(61) 



Here and hereafter we recover the reduced Planck scale Mq. The perturbed equations of 
motion in the longitudinal gauge are given by 



<£> • + HQ* 



Q2y" 

5 Vl1 + 3H5cp t1 + — + — 2 



dV 



= 4^$, - 2-=$, + ^ 



d 2 V 



:5ipij, 



for arbitrary i, j. There is another convenient equation given by 



(62) 
(63) 



*< = E 



-3#0 



A- 2 -2 



(64) 



Before giving the initial conditions for Sipij and $j, we determine the normalization of 
Sifiij. The conjugate momentum of Xi is given by a 3 Xj, which leads to the equal-time commu- 
tation relation [Xi(x, t),Xj(x', £)] = ia _3 (t)5jj5(a; — cc'). Thus the Wronskian normalization 
conditions for Stpij are given by 



(65) 



Then, we find the WKB solutions of Eq. ( |63"D in the short wavelength approximation (k/a ^> 
H), 



-T^exp (-ft/*) 
V2ka V J a J 



(66) 



which give the initial conditions for 5(py. Differentiating these with respect to the cosmic 
time, we obtain the initial conditions for Sip^, 



8cp 



ij — Sij 



<k 1 
M/ 2^ 



2ka 2 



exp I —ik I — 



(67) 
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Note that the exponent can be set to zero because the origin of the conformal time = / dt/a 
is arbitrary. The initial conditions for are derived from Eq. ( p4|) by the use of those of 
Sifij and Sifij. 

Since the vacuum expectation value of the squared of the operator $>a( x , t) is given by 
(o\& A (x,t)\o) = / j0j- 2 J2\Mk,t)\ 2 , (68) 

we define the primordial spectrum /c 3 / 2 $(fc) asf] 

fc 3 /V(0|$ 2 (fe)|0)/5(0) = 2 ( 69 ) 

All quantities are normalized by the combination of M G and if = v/y3. Concretely, 
& = h/M G , V = V/(M G H 2 ), t = tH , fop^ = 5^Hl'\ = $,M G # 1/2 , and k = kH^ 1 . 
Thus all terms except k 2 /a 2 in the equations of motion become of the order of unity, which 
is essentially important for numerical calculations, avoiding rounded errors. Also, in order 
to confirm the results of our numerical calculations we compare the spectrum derived from 
the evolution equation (|62| ) and that obtained from the constrained equation (|64|). Both 
spectra coincide up to the order of 1CT 8 . 

The spectra of primordial fluctuations for the cases of the MACHO PBHs and the BESS 
(GRBs) PBHs are depicted in Figs. p] and ^|. As is easily seen, the numerical results 
reproduce excellently analytic estimates, which is confirmed to be correct. 

IV. DISCUSSION AND CONCLUSIONS 

In this paper we have minutely investigated a natural double inflation model in SUGRA. 
By virtue of the shift symmetry, chaotic inflation can take place, during which the initial 
value of new inflation is set. The initial value of new inflation is adequately far from the local 
maximum of the potential so that primordial fluctuations within the present horizon scale 
are attributed to both inflations. That is, fluctuations responsible for the anisotropy of the 
CMB and the large scale structure are produced during chaotic inflation, while fluctuations 
on smaller scale are produced during new inflation. Due to the peculiar nature of new 
inflation, fluctuations on smaller scale are as large as of the order of unity, which may lead 
to PBHs formation. As examples we consider PBHs responsible for MACHOs and antiproton 
flux observed by the BESS experiment or short gamma ray bursts. We find that if we take 
reasonable values of parameters, such PBHs are produced in our double inflation model. To 
make sure, we also perform numerical calculations and confirm analytic estimates definitely. 



5 Note that the relation between k 3 / 2 &(k) and $>a(N) used in the previous sections is given by 
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TABLES 
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TABLE I. The charges of various supermultiplets of U(l)^ x Z 2 . 
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FIGURES 




FIG. 1. The spectrum of primordial fluctuations which produces the MACHO PBHs is depicted. 



20 



.mordial Po" 



"e: 



SDectrurn 



io- 8 - 



10- 9 - 



-1 11 11 1 nun 1 nun 11 1111 1 11 1 1 11 1 1 111 
- 




111 1 mini 1111111 1 1 1 inn 1 1 11 1111111 1 111 


— 
: 




II 1 1 1 1 III lllll 

< 






□ □ 


E 
- 




□ = 

□ : 






□ □ 

□ — 


- 




D □ \ 
□ □ -_ 






□ = 
□ Z 






□ _ 


= 




□ = 

□ z 
□ - 






n" 


1 1 


1 


1 1 



1 W 1?0 - 1 1 1 : 1 O 2 ! 0=1 01 *L 0^ 0^ O 8 ! Of ¥0 ¥0 ¥0 TO ¥0 f>0 ¥0 ¥0 ¥0 TOTOW 2 

k [h Mpc" 1 ] 



FIG. 2. The spectrum of primordial fluctuations which produces the BESS (GRBs) PBHs is 
depicted. 
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